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ABSTRACT 

We report the detection, with the CANGAROO-III imaging atmospheric Cherenkov telescope array, 
of a very high energy gamma-ray signal from the unidentified gamma-ray source HESS J1614— 518, 
which was discovered in the H.E.S.S. Galactic plane survey. Diffuse gamma-ray emission was detected 
above 760 GeV at the 8.9 a level during an effective exposure of 54 hr from 2008 May to August. 
The spectrum can be represented by a power-law: (8.2 ± 2.2s(at ± ^.bgys) x 10~^^ x {E/1 TeV)"''' 
cm^^ s^^ TeV^^ with a photon index 7 of 2.4 ± Q.3stat ± O-'^sys, which is compatible with that 
of the H.E.S.S. observations. By combining our result with multi-wavelength data, we discuss the 
possible counterparts for HESS J1614— 518 and consider radiation mechanisms based on hadronic 
and leptonic processes for a supernova remnant, stellar winds from massive stars, and a pulsar wind 
nebula. Although a leptonic origin from a pulsar wind nebula driven by an unknown pulsar remains 
possible, hadronic-origin emission from an unknown supernova remnant is preferred. 
Subject headings: gamma rays: observations — ISM: individual(HESS J1614-518) 



1. INTRODUCTION 
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Recent progress with Imaging Air Cherenkov Tele- 
scopes (lACTs) is enabling the exploration of sites 
of cosmic-ray acceleration in our galaxy. Very high 
energy (VHE) gamma rays are produced by the decay 
of neutral pious which arise from interactions between 
the accelerated protons and interstellar matter, or by 
Inverse Compton (IC) scattering and Bremsstrahlung of 
high-energy electrons. For example, VHE gamma rays 
have been detected from young super nova remnants 
(SNRs) such as RX J17 13.7-3 946 (iMuraishi et al. | 
2000; lEnomoto efaLl l2002al: lA haronia n et al.l 
2006 Jl, RX J 0852.0 - 4622 (Katagiri ct all 120051 : 
lAharonian et al.l l2005cl: lEnoinoto et al. 2006bh . and 
RCW86 (jWatanabe et aP 120031 : lAharonian et all 120091 , 
which show p ossible evidence of co sn iic-ray accel- 
eratio n Ce.g.. |Enomoto et al.l l2002al: iMalkov et al.l 
I2005"; 'Uchiyama et al.l |2007t iTanaka et al.l 120081: 
Yamazaki ct al. 2009). In addition, detections of 
VHE gamma rays from pulsar wind nebu lae (PWNe) 
such as the Crab nebula ( Weekes et al.lll989D and Vela X 
nebula (|Aharonian et al.l l2006bt lEnomoto et al.l l2006aD 
have shown that PWNe also play an important role 
in particle acceleration in the Galaxy. Recently, VHE 
gamma-ray emission related to massive stars such as the 
X-ray binary systems PS R B1259— 63 (lAharonian et al.l 
I2OO5I and LS I +61 303 (I Albert et al.ll2006ll. the voung 
open stellar c lusters Cyg 0B 2 (lAha ronian et al"]|2002D . 
Westerlund 2 (lAharonian et a l. 2007), and Westerlund 1 
(|Ohm et al.l 120101 ) have also been reported. Moreover, 
the Galactic plane survey performed by t he H.E.S.S. 
observatory (jAharonian et al.l l2005bL l2006a| ) discovered 
seventeen unidentified VHE gamma-ray sources, includ- 
ing HESS J1614— 518. Today, unidentified sources are 
the largest class of the 123 discovered VHE gamma-ray 
sources, most of which are located in the Galactic 
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lane Ce.g.. lAharonian et alj (|2005dD : lAbdo et"al] 
200l: lAharonian et alj (I2008D: the TeVCAT catalog, 



iTeVCAT ver.3.400l (|20l'll ). is a useful up-to-date, on-line 
resource). In general, the lack of non-thermal electro- 
magnetic radiation from the radio to the X-ray bands 
may be evidence of hadron acceleration because the 
IC scenario requires a lower magnetic field than the 
typical interstellar magnetic field intensity of a few /xG. 
Revealing the possible radiation mechanism(s) of each 
unidentified source is therefore important for identifying 
the origin(s) of cosmic rays. 

H.E.S.S. reported that HESS J1614-518 had a high 
flux level, 25% of the Crab nebula, above 200 GeV with 
a photon index of 2.4 and an elliptical morphology with 
a semi-major axis of 1 4 zfc 1 arcmin and a sem i-minor 
axis of 9 ± 1 arcmin (jAharonian et al.j I2006af) . The 
peak position has an offset of 8 .7 arcmin to the n orth- 
east from th e central position. iLandi et al.l (|2007i ) and 
iRowell et all ([2008) pointed out that HESS J1614-518 
may be associated with the 40 Myr-o ld young open star 
cluster Pismis 22 (jPiatti et al.l |2000[) which is located 
within the VHE gamma-ray emission region at a dis- 
tance of 1.0±0.4kpc and has sufficient luminosity to pro- 
duce the observed gamma-ray luminosity, assuming 20% 
energy conversion from the stellar winds of ten B-type 
stars. However, there are several issues in identifying 
HESS J1614-518 with Pismis 22 since the size of Pis- 
mis 22, 2.0 arcmin in diameter, is one order of magnitude 
smaller than the VHE gamma-ray emission size and the 
location has a 12 arcmin offset from the VHE gamma-ray 
emission peak. In addition, there has been no detailed 
discussion of the radiation mechanism. 

The X-ray satellite Suzaku observed this region 
with the X-ray imaging spectrometer (XIS) and found 
three X-ray sources in their follow-up observation in 
2006 (jMatsumoto et al.l I2008D . One of these, called 
Suzaku source A, is located very close to the VHE 
gamma-ray peak position with an offset of 0.8 arcmin. 
The spectrum is well fitted by a single power-law 
model with a photon index of I.TSIq'jo ^^"^ ^ hydrogen 
equivalent column density of 1.21^q4° x 10^^ cm~^. 
The distance to Suzaku source A is approximately 
10 kpc, which was derived from the hydrogen equivalent 
column density using the total Galactic H i column 
density toward HESS J161 4-518 of - 2.2 x 10^2 cm^^ 
(jDickev fc LockmanI |1990[) . The size of the X-ray 
emission region is slightly larger than the Suzaku 
Point Spread Function (PSF) of 1.8 arcmin and smaller 
than the size of the VHE gamma-ray region. This 
difference is seen in PWNe such as HESS J182 5-137 
(lAharonian et al.l l2006dl: lUchiyama et al.l I2009D and 



Vela X (jMarkwardt fc Qgelmanl 119951: lAharonian et al.l 
l2006bD and could be explained by the difference between 
the synchrotron cooling time of the electrons that 
radiate X-rays and those that produce TeV gamma rays. 
Electrons with an energy of 100 TeV radiate X-rays and 
immediately lose their energy by synchrotron cooling 
(e.g., the energy- loss timescale is ~ 10^ yr ass uming a 
magnetic field of 20 /xG (jSturner et al.ril997D ). while 
electrons with an energy of 1 TeV, which are responsible 
for the VHE gamma ray emission through IC scattering, 
are more slowly cooled by synchrotron radiation (e.g., 
the energy-loss timescale is '^ IC* yr assuming the 



same parameters as above) and can travel further from 
their source. However, the ratio between the observed 
VHE gamma-ray and X-ray fluxes, F(l-10 TeV) / 
F(2-10 keV) of ^34, is much larger than those of 
known PWNe — 2.6 x IQ-^, 0.7, and 1. 5 for the Crab 



MSH 1 5-52, and Vela X resp e ctivelv (fCaensler et al.l 
19991: IWiUingale et all l2nnit IGaensler et al.l 12002 : 



Dodsqn et al.l 120031: lAharonian et al.' '200i |2005a, 
2006bl: iManzah et al.l [2007: Nakamori et al. 2 00l). Nev 



erthe less, recent studies of HESS J1640-465 (Funk et al.l 
I2OOI and HESS J1804-216 (Higashi et al. 2008) claim 
that this large ratio can be explained by a time-evolving 
electron injection model, in which the number of 
electrons injected into space by the pulsar decreases 
proportionally to the spin-down of the pulsar. On the 
other hand, this large ratio is also expected in an old 
SNR with an age of ^ 10^ yr, because of the difference 
between the cooling times of electrons and protons 
(jYamazaki et al.l 120061 ). We therefore discuss both a 
PWN scenario and an SNR scenario in this paper. 

Suzaku source B is positioned towards the center of 
HESS J1614— 518 and is coincident with the position of 
Pismis 22. Since the hydrogen equivalent column density 
derived from the Suzaku spectrum is (1.1 ± 0.21) x 10^^ 
cm~2^ which is comparable with that of Suzaku source A, 
Suzaku source B may lie at a similar distance as for 
Suzaku source A. This source has a non-thermal X-ray 
emission with a photon index of 3.19 ± 0.32. This soft 
index and X-ray luminosities of 7.7 x 10'^'* ergs s~* and 
4.5 X 10^^ ergs s'^ in the 2-10 keV and 0.5-10 keV 
ranges, respectively, assuming a distance of 10 kpc, are 
typic al values for a n anomalous X-ray pulsar (AXP) 
(Fahl man fc Gregor v 1981; Kuipcr ct al'. 2006). The pos- 
sible existence of the AXP also suggests this source as 
an SNR, since AXPs are usually associated with SNRs, 
e.g., IE 2259-^586 with CTB 109 (Fahlman & G regorv 



1981) and IE 1841-045 with Kes 73 ( Vasisht fc GottheU 
1997() . Since the position of this source is coincident 



with Pismis 22, there is a possibility that this emission 
originates from the stellar winds from the stellar clus- 
ter. Non-thermal X-ray emissio n from a stellar c luster 
was reported from Westerlund 1 (jMuno et al.ll2006l ) , and 
TeV gamma rays were recently detected from this object 
(jOhm et al.ll2009l ). However, this positional correlation 
may be only a chance coincidence since the estimated 
distances to Suzaku source B and Pismis 22 are different 
by an order of magnitude. Although Suzaku source B 
might be marginally extended, it is diflricult to quanti- 
tatively estimate the spatial extension with the Suzaku 
PSF of 1.8 arcmin. If Suzaku source B is actually ex- 
tended, additional scenarios besides an SNR could be 
co nsidered, e.g., a PWN fro m a pulsar/AXP as discussed 
in iMatsumoto et aLl (|2008l ) , or emission from the unre- 
solved hot stars in Pismis 22. 

The other source , Suzaku source C, is a late B-type 
star as described in IMatsumoto et al.l (|2008l) , and thus is 
not a possible counterpart of HESS J1614— 518. 

Swift observed this region with the X-ray telescope 
(XRT) and foun d six X-ray sourc e s (he reafter Swift 
sources 1 to 6) (|Landi et al.l 120061 I2q07| ). All these 
sources were point-like and no diffuse emission was found. 
Two sources. Swift sources 1 and 4, are located within 
the field of view (FOV) of the Suzaku observation. Swift 
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source 1 is located close to Pismis 22 with an off- 
set of 42 arcsec. This source is also coincident with 
Suzaku source B. Swift source 4 is coincident with Suzaku 
source C. Swift sources 1, 2, 3, 5 are probably stars, while 
the nature of Swift sources 4 and 6 were not identified, 
probably due to the poor statistics. Although Suzaku 
source A was located in the FOV of the Swift XRT, it 
was not detected with Swift probably due to the limited 
exposure time (~1700 s) and/or the small effective a rea. 
The Fermi-LAT collaboration (|Abdo et al.ll2010al ) re- 
ported the detection of gamma rays in the 100 MeV 
to 100 GeV band from IFGL J1614.7-5138c positioned 
2.7arcmin away from the peak position of the VHE 
gamma-ray emission. In the radio band, no counter- 
part has been found in the HESS J1614— 518 region; 
there is no enhancement in the 843 M Hz band, where 
the rms noise level is ~2 mJy arcmin^^ ( Bock et al.lll999l : 
iMurphv et al.l [2007|). In this paper, we present TeV 
gamma-ray observations of HESS J1614— 518 with the 
CANGAROO-III telescopes and discuss the possible 
counterpart and the radiation mechanism by consider- 
ing multi-wavelength observations. 

2. OBSERVATIONS 

CANGAROO-III is an array of four lACTs (Tl, T2, 
T3, and T4), located at Woo mera, South Australi a 
(136°47E, 31°06S, 160 m a.s.L) (lEnomoto et al.ll2002bn . 
The oldest telescope, Tl, which was the CANGAROO- 
II telescope, has not been in use since 2004 due to 
its smaller FOV and higher energy threshold. Each 
telescope has a 10 m diameter reflector which con- 
sists of 114 segmented fiber reinforced p lastic spherical 
mirro rs mounted on a parabolic frame (jKawachi et al.l 
1200 It ). The imaging camera system consists of 427 
photomultipliers (P MTs) and has a FOV of 4.0 deg 
(|Kabuki et alll2003f ). The PMT signals are digitized by 
charge analog-to-digital converters ( ADCs) and mult i- 
hit time-to-digital converters (TDCs) (jKubo et al.ll2001[ ). 
The observations were carried out from 2008 May to Au- 
gust using a wobble mode in which the pointing position 
was shifted both in declination and right ascension be- 
tween ±0.5 deg from th e target position every 20 minutes 
(jNakamori et al.ll2008() . The target position was (R.A., 
decl. [J2000])=(243°.579,-51°.820) which is the center 
of the source position reported by H.E.S.S. 

The light-collecting efficiency, including the reflectiv- 
ity of the mirror segments and the light guides, and the 
quantum efficiency of PMTs, was monitored by a muon- 
ring analysis with individual trigge r data taken in the 
same period ([Enomoto et al.|[2006aD . The average quan- 
tity of light per unit arc-length of muon rings is approxi- 
mately proportional to the light-collecting efficiency. We 
also did not use the second oldest telescope, T2, since 
it was very difficult to calibrate with muon-ring data to 
calculate the efficiency at t his time due to the d eterio- 
ration of mirror reflectivity (jEnomoto et al.ir2009f) . The 
two telescopes T3 and T4 were used. From the muon- 
ring analysis for the data taken in this period, the light- 
collecting efficiency of each telescope, which is used in 
the Monte Carlo simulations, with respect to the orig- 
inal mirror production time, was 0.58 and 0.50 for T3 
and T4, respectively. We reject the data for which the 
average trigger-rate over a one minute period was under 
5 Hz to remove data taken in cloudy conditions. The ef- 



fective exposure time amounts to 53.6 hr and the energy 
threshold was 760 GeV. 

3. ANALYSIS 

The st andard analysis of the CANGAROO-III co llab- 
oration (jEnomoto et all l2006at iKabuki et al.ll2007[) was 
applied to the data. The calibrations for the cameras 
and ADCs were carried out daily using LEDs. After cal- 
ibration, the recorded charges of each pixel in the camera 
were converted to the number of photo-electrons. Time- 
walk corrections for TDC data were carried out using 
data taken by changing the luminosity of the LEDs. Af- 
ter calibration, every shower image was cleaned through 
the following criteria. Only pixels that had >5.0 photo- 
electrons were used as "hit pixels" . Clusters of five or 
more adjacent hit pixels with arrival times within 30 nsec 
of the average hit time of all pixels were recognized as 
a shower cluster. Before calc ulating image moments — 
the "Hillas parameters" (iHill as 1985 ) — we ap plied an 
"edge cut" to the data (jEnomoto et al.|[2006b[ ) and re- 
jected events having hit pixels in the outer-most layer of 
the camera. The orientation angles were determined by 
minimizing the sum of the squared widths with a con- 
straint given by the distance predicted by Monte Carlo 
simulations. 

In order to derive the gamma-ray likeliness, we 
used the Fisher Dis criminant method (Fis her I 119361 : 
lEnomoto et al.l[2006a| ). The input parameters were 

P={W3,Wi,L3,Li), 

where W^.W^.L^.L^ are energy-corrected Widths and 
Lengths for T3 and T4 camera images. The Fisher Dis- 
criminant (hereafter FD) is defined as FD = a-P, where 
Q? is a set of coefficients mathematically determined in 
order to maximize the separation between two FD dis- 
tributions for gamma rays and hadrons. 

To calculate the background, we selected a ring re- 
gion around the target, 0.3 < 9^ < 0.5 deg^, where 9 is 
the angular distanc e to the center of HESS J1614— 518 
(jHigashi et al.ll2008[ ). We obtained the FD distributions 
of hadrons Ff, from this region and gamma rays Fg from 
Monte Carlo simulations. Finally, we were able to fit the 
FD distribution of the events from the target with a lin- 
ear combination of these two components. The observed 
FD distributions F were represented by 

F = pFg + {l- p)Fb, 

where /3 is the ratio of gamma-ray events to the total 
number of events. Here /3 is a fitting parameter and the 
obtained FD distributions are shown in Figure [T] 

4. RESULTS 

The obtained 9'^ plot is shown in Figurej^jwith the PSF 
of our telescopes. Above 760 GeV we detected 950±107 
excess events within 9^ < 0.2 deg^. The morphol- 
ogy of gamma-ray-like events, obtained from a gaussian 
smoothing with the CANGAROO-III PSF of 0.24 deg, 
is shown in Figure IH The extent of the VHE gamma- 
ray emission was estimated by a 2D Gaussian fit on 
our excess map. The obtained standard deviation was 
0.44 ± 0.03 deg which is broader than the CANGAROO- 
III PSF. The centroid position was determined to be 
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Fig. 1. — FD distribution. The black circles show the FD ob- 
tained from the ON source region, 9^ < 0.2 deg^. The red and blue 
lines are the background and gamma-ray component estimated by 
the fit procedure described in the text. The green circles are ob- 
tained by the subtraction of the background from the ON source 
region. 

(R.A., decl. [J2000])= (243°. 634, -51°.950). The offset 
from the best-fit position reported by H.E.S.S. is (AR.A., 
Adecl.) = (0°.055 ± 0°.018 -0°.130 ± 0°.033). The off- 
set is not significant compared with our PSF. A system- 
atic difference due to the difference in energy thresholds 
between H.E.S.S. and CANGAROO-III may also con- 
tribute to this offset. The reconstructed VHE gamma- 
ray differential flux is shown in Figure |4] together with 
the H.E.S.S. measurements. The spectrum is compat- 
ible with a single power-law: (8.2 ± 2.2stat ± 2.5sys) x 
10-12 X (E/1 TeV)-''' cm-2 s'^ TeV'^ with a photon 
index 7 of 2.4 ± 0.3stat ± 0.2sys- The relevant systematic 
errors are due to the atmospheric transparency, night- 
sky background fluctuations, uniformity of camera pixels, 
and light-collecting efficiencies. In addition, to estimate 
the systematic error due to the size of signal integra- 
tion region, we changed the region from 9^ < 0.14 deg^ 
to 0.30 deg^, which was included in the systematic er- 
ror. For comparison with the ring-region background, we 
took background events from opposite positions of HESS 
J1614— 518 observations in the wobble mode, and then 
obtained a differential fiux of (6.4 ± 2.Qstat ± 2.4:sys) x 
10^12 X (E/l TcV)"'^ cm-2 s~^ TeV^^ with a photon 
index 7 of 2.4 ± O.Qstat i 0.3sjys, which was consistent 
with that derived with the ring-region background. The 
VHE gamma-ray extension, centroid position, and flux 
obtained with CANGAROO-III are consistent with re- 
sults from H.E.S.S. This result suggests that the VHE 
gamma-ray emission was unchanged between 2004 and 
2008. 

5. DISCUSSION 

We now discuss the plausible radiation mechanisms 
of HESS J1614-518 using the resuhs of CANGAROO- 
III, H.E.S.S., Fermi, and Suzaku observations. Since 
the spectra of Swift sources were not available, we did 
not use the Swift data. Figure [5] shows the morpho- 
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Fig. 2. — The 8^ plot, where 9^ = corresponds to the fit ted cen- 
ter of gravity of HESS J1614-518 from H.E.S.S. (Aharonia n et al.l 
I2006al ). The blue data points represents the excess events in each 
9^ bin and the red solid line represents our PSF derived from the 
Monte-Carlo simulation. 
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Fig. 3. — Morphology of gamma-ray-like events. The number of 
excess events per 0°.04 X 0°.04 cell is smoothed by a Gaussian with 
cr = 0.24 degree, which is the CANGAROO-III PSF, and plotted 
in equatorial coordinates. The black solid contours show the VHE 
gamma-ray emission seen by H.E.S.S. Lines correspond to 20, 30, 
40, 50, & 60 gamma-ray counts. The red cross shows the H.E.S.S. 
center of gravity of HESS J1614-518 (Aharonian et al. 2006at) . 

logical relationship between each observation. The non- 
thermal X-ray emission from Suzaku source A is posi- 
tioned very close to the H.E.S.S. gamma-ray peak, the 
position of the Fermi source IFGL J1614.7-5138c, and 
within the emission region detected with CANGAROO- 
III. Thus, this could be the most likely counterpart for 
HESS J1614-518. We note here that since the FOV 
of the Suzaku observation covered only the part of the 
TeV gamma-ray emission region, as shown in Figure [SJ 
the current observed X-ray flux may only be a fraction 
of the entire X-ray emission from the entire region of 
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Fig. 4.— Differential flux of HESS J1614-518. Squares and 
circles show the CANGAROO-III and the H.E.S.S. data points, 
respectively. The best fit power-law from this work is shown by 
the dotted line. 

the VHE gamma-ray emission. To discuss the emission 
mechanism with more accuracy, further X-ray observa- 
tions of the entire region of VHE gamma-ray emission are 
needed. For further constraints on emission models, we 
tried to derive the flux at 8 ^m at the Suzaku source A 
position, from t he archival data w ith the Infrared Array 
Camera (IRAC. [Fazio et all ()2004[) 1 onboard the Spitzer 
Space Telescope. However, since there was contamina- 
tion from a nearby source, we obtained only an upper 
limit. 

SNR Scenario: Suzaku source B, which may be related 
to a possible AXP, is positioned roughly in the center 
of the VHE gamma-ray emission. We thus postulate a 
scenario in which a supernova explosion occurred at the 
position of Suzaku source B and the shock of the SNR has 
now reached the position of Suzaku source A, emitting 
both the X-ray and gamma-ray emission. 
PWN Scenario: A PWN could also emit dif- 
fuse gamma-ray emission. Five pulsars have 
been found in this region, PSR J1611— 5209, 
PSR J1612- 

and PSR 

lATNF Pulsar Catalogue ver. 1.381 ). As described m 
section 1, Suzaku source B would be a PWN if associ- 
ated with a pulsar or AXP. The smaller size of Suzaku 
source A or source B compared to that of the VHE 
gamma-rays also appears in other PWNs because of 
synchrotron cooling (section 1). 

We will discuss the PWN scenario and the associated 
pulsar which could supply enough particles to reproduce 
both the X-ray emission of Suzaku source A or source B 
and the VHE gamma-ray emission. 
Stellar Wind Scenario: The young open cluster Pis- 
mis 22 is located towards the center of HESS J1614— 518 
and is also a possible counterpart. Its age is '^ 4.0 x 10^ yr 
and the distance is l±0.4kpc from the Earth. The coinci- 
dence between a young open cluster and a VHE gamma- 
ray source is also seen in Westerlund 2 and Cyg OB 2 as 
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Fig. 5. — Morphological relationship between the X-ray and 
gamma-ray observations. Blue, magenta, and green contours 
show the emission regions of CANGAROO-III, H.E.S.S., and 
Suzaku, respectively. Orange dashed circles show the 68% and 
90% error ellipses of the position of the Fermi source. The 
red circle shows the Pismis 22 position and purple squares 
show the positions of nearby p ulsars ((Manchester et al. 2003 
lATNF Pulsar Catalogue ver. 1.381) . The bold bla ck line shows the 
observed region by the Suzaku XIS in Matsumot o et al.l II2008I ) ex- 
cluding the calibration source region on the corners. 

described in section 1. Stellar winds from massive stars 
could form a shock front, accelerate c harged particles, 
and produce the high energy radiation ()Voelk fc FormanI 
[Mi [Bedn"arek"200r). 

Another possibility is an association between an SNR, 
PWN and the open cluster, since Pismis 22 is old enough 
for some massive stars to finish their life as supernovae. 
In addition, as described in section 1, binary systems 
that emit VHE gamma rays also have been discovered. 
However, HESS J1614— 518 does not seem to be associ- 
ated with a binary system, since all sources of this type 
have point-like emission, which is in contrast to the re- 
sults of CANGAROO-111 and H.E.S.S. We discuss the 
above scenarios in detail in the following subsections. 

5.1. SNR Scenario 

In the SNR scenario, we assume that the X-rays from 
Suzaku source A and the VHE gamma rays are emitted 
by charged particles accelerated by the shock in the SNR 
shell, and Suzaku source B is an associated AXP which is 
positioned in the center of the SNR shell. Additionally, 
we will discuss the possible correlation between the SNR 
and Pismis 22. 

First we examine a leptonic model to explain the ob- 
served SED (Fig. [S]). For the X-ray spectrum, we 
use the Suzaku spectrum of Suzaku source A corre- 
lated with the VHE gamma-ray peak, with a statisti- 
cal error at the 90% confidence level. For the Fermi 
spectrum, we used a 2a statistical error and system- 
atic er rors of 1.8g in the flux and 1.2a in the photon 
index (jAbdo et al.l l2010al) . Figure [Hj shows the simple 
one-zone leptonic model curves for HESS J1614— 518. 
Here we calculated synchrotron, inverse Compton, and 
Bremsstrahlung model curves for a single power-law with 
an exponential cutoff electron spectrum, dNe/dE,, — 
Ki,E~^'' exp{—Ee/Emax_e), whcrc Kf, is the normaliza- 
tion factor, Ee is the electron energy, Fe is the spec- 
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tral index of the electrons, and Emax^e is the maxi- 
mum accelerated electron energy. To calculate the IC 
radiation, we used cosmic-ray microwave background 
and an interstellar radiation fiel d (ISRF) derived from 
the GALPROP package (v 50p) (jPorter fc StrondlMia 
IStrong &: Moskalenkoll2006|) to estimate the seed photon 
field around the HESS J1614— 518 region. Energy densi- 
ties of 1.4 eV cm~'^ and 5.5 eV cm""^ for IR and optical 
light were obtained, respectively. In addition, this radi- 
ation field changed by less than an order of magnitude 
when we varied the distance from Ikpc to lOkpc, with 
values in the range from 0.9 to 1.7 eV cm~^ and 1.1 to 5.5 
eV cm~^ for IR and and optical light, respectively. We 
fixed the power-law index to Fe = 2.0 and fitted the VHE 
gamma-ray spectrum by the IC emission. The maximum 
energy and total energy of the electrons obtained were 
4.2±1.5 TeV and 1.9 x lO''^ x {d/Wkpcf ergs. From com- 
parison between the synchrotron model and the Spitzer 
upper limit, an upper limit of a magnetic field is de- 
termined to be 6 /iG. The harder and fainter spectrum 
in the X-ray band compared to the VHE gamma-ray 
spectrum did not allow the observed X-ray and VHE 
gamma-ray spectra to be produced with synchrotron 
and IC emission, respectively, by a single power law 
distribution of electrons. There is also the possibility 
that Bremsstrahlung produces both the X-ra y and VHE 
gamma-ray emission (' Uchivama et al.l l2002f) , as shown 
in Figure [SI This model gives a good reproduction for 
an ambient m atter density Up of 600 p cm~^. However, 
IRowell et al.l (|2b08) reported that they found no obvious 
overlapping molecular clouds across a range of inferred 
distances up to ~6 kpc, with the NAN TEN i2cO(J=l- 
0) survey data (jMatsunaga et al.ll200lD . In addition, we 
estimated an ambient matter density from the velocity- 
integrated data of the CO survey (Dame et al. 20_0_1.) to 
be ~80 or 8 p cm~'^ for 10 kpc and 1 kpc, respectively. 
Thus, the Bremsstrahlung model that requires an am- 
bient matter density of 600 p cm~^ was rejected. Fur- 
thermore, there is difficulty in explaining the morpholog- 
ical difference between the X-ray and the VHE gamma- 
ray emission because it requires an unlikely situation in 
which the relatively high energy (multi-TeV) electrons 
that are responsible for the VHE gamma-ray emission 
are distributed over a more extended region than the 
relatively low energy (multi-keV) electrons which are re- 
sponsible for the X-ray emission. 

We also checked that our estimate of the non-thermal 
synchrotron flux did not violate the thermal optical emis- 
sion. The intrinsic optical background flux over the solid 
angle for the gamma-ray emission region ^ 5 x 10^ eV 
cm~^s~^ is obtained from Uopt{c/4:TT)fl, where Uopt is the 
energy density of 5.5 eV cm~^ for the optical light, c is 
the speed of the light, and fl ^ 7r(0.2 x 7r/180)^ sr is the 
solid angle for the region. The optical background radia- 
tion is 10^ times larger than the optical synchrotron flux 
for a magnetic field of 6 ^G (Fig. ^. 

Second, we examined a neu tral-pion decay model. 
Based on a model proposed bv lYamazaki et all (|2006[ ). 
only nucleonic particles remain in an old SNR with 
an age of ~ 10^ yr, while primary electrons have al- 
ready lost most of their energy by the synchrotron cool- 
ing. Figure [7] shows the SED with the assumption 
that the population of accelerated protons can be ex- 



pressed by a single power-law with an exponential cut- 

off, dNp/dEp — KpEp " exp{—Ep/Emax.p)- We set the 
power-law index to Tp = 2.0. The best-fit cutoff en- 



ergy was obtained to be E„ 



= 36 ± 18 TeV. The 



total energy of high energy protons was calculated to be 
Wp = 1.2 X 1052(np/l pcm-3)-i(d/10kpc)2 ergs. By 
setting Up = 100 p cm~'^, the efficiency of energy con- 
version to accelerate protons is 10% for a typical to- 
tal supernova explosion kinetic energy of ^ 10^^ ergs. 
As described above, no obvious molecular cloud was 
found in the NANTEN 12CO(J=1-0) survey data. In 
addition, an ambient matter density from the velocity- 
integrated data of the CO survey (Dame et al., 2000) 
was estimated to be ~80 or 8 p cm~'^ for 10 kpc and 1 
kpc, respectively. Thus, an assumption of 100 p cm^'^ 
is likely for 10 kpc, but not for 1 kpc. Further ob- 
servations are necessary to investigate the validity for 
such assumed density. Assuming the spectral index, 
Fp = Fe = 2.0, the maximum energy of primary elec- 
trons is determined to be Emax^e £ Emax.p- The 
turnover energy of synchrotron emission is determined 
from Eturn = 22 keV x(£;„a._e/50 TeV)2 x {B/200fiG). 
Since the hard index of the X-ray spectrum required 
Eturn > 10 keV, a lower limit for the magnetic field of 
B >200 /iC was determined from this equation. The 
model curve for this lower limit condition was shown in 
Figure [T] From this magnetic field value, we set a lower 
limit to the number ratio of protons to primary electrons, 
Kpe = Kp/Ke > 2.1 X I0^{np/1 p cm-3)-i. 

We also calculated the contribution of emissions from 
secondary electrons from p-p interactions between the 
same proton population as above and the ambient mat- 
te r density of lOOpcm ^ '^. W e followed the calculation 
in iKelner fc AharonianI ()2008D to derive the spectrum of 
the secondary electrons. Assuming a distance of 10 kpc 
from the Earth, the distance between Suzaku source A 
and source B was calculated to be 35 pc. Thus, we as- 
sumed that the radius of the SNR is 35 pc and roughly 
estimated the a g e of t he SNR using the equation (2) in 
lYamazaki eTan (|2006[ ) to be 3 x 10'' yr. Thus, the emis- 
sions from the secondary electrons were derived by as- 
suming continuous injection of electrons produced by a 
constant proton s pectrum over 3 x 10^ yr with a magnetic 
field of 200 AiG (JAtovan fc AharonianI Il99l . The ob- 
tained synchrotron curve is also shown in Figure [7l The 
obtained inverse Compton and Bremsstrahlung emissions 
were able to be neglected since the number of electrons 
is sufficiently small. Since the synchrotron emission from 
the secondary electrons was not able to explain the X-ray 
emission, the X-ray emission might originate in the syn- 
chrotron emission from the primary electrons, as shown 
in Figure El or other emission mechanisms. Since the 
synchrotron emission from the secondary electrons dom- 
inates below the infrared band, the detection of the emis- 
sion in the radio to infrared bands could support the hy- 
pothesis that the VHE gamma-ray emission is produced 
by the neutral-pion decay. Because the Spitzer upper 
limit was above the predicted flux of synchrotron emis- 
sion, more detailed observations are needed. 

To discuss the possible association between the SNR 
and Pismis 22, we estimated the SNR age using the dis- 
tance to Pismis 22 with the assumption that a super- 
nova explosion occurred at the position of Pismis 22 and 
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Fig. 6.— SED with leptonic model curves for HESS J1614-518. 
The dash-dotted and soUd blue lines show IC and synchrotron emis- 
sion derived from the single power-law electron spectrum with an 
exponential cutoff to fit the VHE emission with a magnetic field of 
6 fiG. The dashed black line shows a Bremsstrahlung curve for a 
number density of ambient matter of 600 p cm~^. 

the shock front has now reached at the VHE gamma- 
ra y peak position . By using the same equation (2) 
in lYamazaki et all (|2006l ). the age was obtained to be 
3 X 10^ yr for a distance of 1 kpc. Since the lower hmit 
of a magnetic field was obtained to be B >200 /iG as 
above, the synchrotron cooling time of 100 TeV electrons 
decreased to 1 yr. For a distance of 1 kpc, the required 
total energy of protons Wp could be reduced if the num- 
ber density of ambient matter Up is the same as for a dis- 
tance of 10 kpc. If the total energy of protons Wp is fixed 
to be 10^° ergs, the required density of ambient matter 
is reduced to Up — 1 p cm^'^ for a distance of 1 kpc. 
This value was comparable with the typical number den- 
sity in the interstellar field and does not contradict the 
fact that no obvious molecular cloud was found in the 
NANTEN 12CO(J=1-0) survey data. The contribution 
of emissions from secondary electrons was shown in Fig- 
ure [71 assuming the injection time of 3 x 10^ yr with a 
magnetic field of 200 /iG. The synchrotron emission from 
the secondary electrons was not able to explain the X-ray 
spectrum of Suzaku source A. The X-ray emission might 
originate in the synchrotron emission from the primary 
electrons, or other emission mechanisms. Since the flux 
of synchrotron emission from the secondary electrons in 
the radio to infrared bands was lower than that of the 
case of a distance of 10 kpc, as shown in Figure [71 a de- 
termination of the spectrum below the infrared band is a 
key to reveal the origin of the SNR. Additionally, an ob- 
servation of thin thermal plasma in the X-ray band will 
provide useful information such as plasma temperature 
or chemical abundances. For example, the detection of 
high abundance of a-elements (O, Mg, Si, S, Ca, and Ti) 
compared to th at of iron, which is exp ected in a massive 
star explosion ([Kobavashi et al.l 120061. may support the 
SNR scenario. In fact, a recent Suzaku observation of 
the open cluster Westerlund 2 detected metal-rich ther- 
mal emission, suggesting that the diffuse X-ray and VHE 
gamma-ray emission may have o riginated from a hyper- 
nova remnant ^Fuiita et al.ll2009l ). 

5.2. Stellar Wind Scenario 

The VHE gamma-ray emission might have been pro- 
duced by hadrons accelerated in winds from massive stars 
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Fig. 7. — SED and the model curve for neutral pion decay (blue 
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in Pismis 22 ("Voclk fc Formanlll982HBedn arck'200r). A 
fraction of the stellar wind energy can be transferred to 
relativistic particles. Assuming that the shock acceler- 
ation generates a single power-law spectrum of primary 
particles, we can apply the discussion made in the SNR 
scenario. We discuss the energetics for the hadronic ori- 
gin here to produce the observed gamma-ray emission. A 
single 0-type star loses mass at a rate of M = 10~^Mq 
per year with a ste llar wind velocity of ~ 1500 km s 
(jCastor et al.lll975l) . The rate of kinetic energy emitted 
from the single star is 7 x 10^^ erg s~^ . If we assume 
an energy conversion efficiency to particle acceleration 
of 5 %, which is t he maximum effi ciency adopted for a 
hadronic model in iBednarekl (|2007[ ). an ambient matter 
density of 100 p cm"'', a distance of Ikpc, and an age 
of 40Myr, two 0-type stars are required in Pismis 22 to 
produce the observed VHE gamma-ray spectrum in the 
hadronic scenario. However, no obvious molecular cloud 
has been found at this distance in the NANTEN data, 
as described above. 

5.3. PWN Scenario 

PWNe are the largest class of identified Galactic 
VHE gamma-ray sources. We discuss the possibility of 
HESS J1614-518 b eing a PWN ii i this subsection. 

As calculated in iRowell et al.l (|2008[ ). the spin-down 
luminosity of each of the five nearby pulsars (Fig. [Sj 
is smaller than the TeV gamma-ray luminosity. Thus, 
none of already known pulsars can be associated with 
HESS J1614-518. There is also the possibility that an 
undiscovered pulsar with a high spin-down power, suf- 
ficient to explain the observed gamma-ray luminosity, 
might be located in the vicinity of Suzaku source A or 
source B. We estimated the pulsar age to be 24 kyr 
and 23 kyr for Suzaku source A and B, respectively, 
using the correlation between pulsar age and the ra- 
tio of gamma-ray flu x and X-ray flux from the PWN 
(jMattana et al.l [20091 ). For this relatively old age, we 
should apply the time-evolving electron injection model 
(section 1). We applied thi s model to HESS J16 14-518, 
following the calculation in iHigashi et al.l (|2008f) applied 
to HESS J1804-216. Although this model can explain the 
large ratio between the X-ray and TeV gamma-ray fluxes 
assuming a single power-law electron distribution, the 
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model showed a very large discrepancy with the sub-GeV 
flux observed with Fermi, the model curve is 30 times 
larger than the observed flux at 0.1 GeV. Thus, the time- 
evolving electron injection model with a single power-law 
electron distribution was rejected. The Me V/ GeV com- 
ponent could arise from the different mechanism than the 
TeV emission. 

Although the present sensitivity in the radio band may 
not be sufficient to detect this unknown pulsar, further 
observations in the GeV band could detect a radio quiet 
pulsar li ke Geminga, which was detected with CGRO 
EGRET fBertsch et al.|[T99^ . or sixteen previously un- 
known pulsars whi ch were recently discovered with Fermi 
(|Abdo et al.ll2009t ). In addition, future X-ray observa- 
tions could detect pulsed emission from Suzaku source A 
or source B. G iven these detections, emission models 
for PWNs (e.g.. lAbdo et al.! '2010b: 'Tanaka fc Takaharal 
120101: ISlane et al.l I2010t IBucciantini ct al. 20TT1 ). using 
a broken power-law distribution of electrons, might be 
able to reveal that the VHE gamma-ray emission of this 
unidentified source originates from a PWN. 

6. CONCLUSION 

The observation of HESS J1614-518 with the 
CANGAROO-III telescopes confirms the VHE gamma- 
ray emission reported by H.E.S.S. The differential energy 
spectrum can be fitted with a single power law: (8.2 ± 
2.2stat ± 2.bsys) X 10-12 X (^/lTeV)-T cm-2s-iTeV-i 
with a photon index 7 of 2. 4±0.3stat ±0.2^2^5. We discuss 
the possible counterparts for this object using the results 
of observations with Suzaku and Fermi. For the SNR sce- 
nario, a one-zone leptonic model was not able to account 
for the observed SED. Hadronic models gave a good re- 
production of the SED and the typical SNR explosion 
energy of ^ 10^^ ergs is able to supply the total energy of 
protons. Since the required number densities of the am- 
bient matter were rip = \ p cm"^ and rip = 100 p cm"'^ 
for a distance to the SNR of 1 kpc and 10 kpc, respec- 
tively, detailed molecular observations could determine 
whether the SNR originated from Pismis 22 (d'-~^l kpc) 
or a farther distance. As there were also differences in 
the spectrum of the emission from the secondary elec- 
trons, a determination of the spectrum below the in- 
frared band would help determine the likelihood of an 
SNR origin. For the PWN scenario, the nearby known 



pulsars are not responsible since the spin-down powers 
are insufficient to produce the observed TeV gamma-ray 
luminosity. Further observations to search a pulsar are 
necessary to investigate the PWN scenario. For the stel- 
lar wind scenario, Pismis 22 was required to contain two 
0-type stars through its entire age from energetics con- 
siderations. However, the required number density of the 
ambient matter of rip = 100 p cm""^ may not be consis- 
tent with the results of the NANTEN observations. 

To identify HESS J1614-518, more detailed multi- 
wavelength observations are required. To discuss the 
stellar wind origin in more detail, a determination of 
the number of OB stars is necessary. For the SNR sce- 
nario, the ultra-high energy resolution of the SXS on- 
board Astro- H could detect line emissions with a high 
abundance of a-elements compared to that of iron, which 
would indicate that HESS J1614-518 is an SNR. To show 
that the Suzaku source B is an AXP, supporting the SNR 
or PWN scenario, a high time-resolution X-ray observa- 
tion is needed to detect a pulsed signal from the source. 
More detailed gamma-ray spect roscopy wit h Fermi, and 
the Cherenkov Telescope Array ()CTAI|2010( ) could deter- 
mine the origin of the accelerated particles. 
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